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Abstract 

We discuss a possible spectral realization of the Riemann zeros based on the Hamiltonian H = xp 
perturbed by a term that depends on two potentials, which are related to the Berry-Keating 
semiclassical constraints. We find perturbatively the potentials whose Jost function is given by 
the zeta function £(<r — it) for a > 1. For a = 1/2 we find the potentials that yield the smooth 
approximation to the zeros. We show that the existence of potentials realizing the zeta function at 
a = 1/2, as a Jost function, would imply that the Riemann zeros are point like spectrum embedded 
in the continuum, resolving in that way the emission/spectral interpretation. 

PACS numbers: 02.10.De, 05.45.Mt, ll.10.Hi 
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I. INTRODUCTION 



One of the most important problems in Mathematics is the proof of the Riemann hy- 
pothesis (RH) which states that the non trivial zeros of the classical zeta function all have 
real part equal to 1/2 [l, 2]. The importance of this conjecture lies in its connection with 
the prime numbers. If the RH is true then the statistical distribution of the primes will be 
constrained in the most favorable way. Accordingto Michael Berry the truth of the RH 
would mean that "there is music in the primes" 3, 41. Otherwise, in the words of Bombieri, 

n 

the failure of the RH would create havoc in the distribution of the prime numbers [5|. 

In so far, the proof of the RH has resisted the attempts of many and most prominent 
mathematicians and physicist for more than a century, which explains in part its popularity 



0, 0, fl. 



For these and other reasons the RH stands as one of the most fundamental 
problems in Mathematics for the XXI century with possible implications in Physics. In fact, 
physical ideas and techniques could probably be essential for a proof of the RH 0, Q . This 
suggestion goes back to Polya and Hilbert which, according to the standard lore, proposed 
that the imaginary part of the non trivial Riemann zeros are the eigenvalues of a self-adjoint 
operator H and hence real numbers. In the language of Quantum Mechanics the operator 
H would be nothing but a Hamiltonian whose spectrum contains the Riemann zeros. 

The Polya-Hilbert conjecture was for a long time regarded as a wild speculation until the 
works of Selberg in the 50's and those of Montgomery in the 70's. Selberg found a remarkable 
duality between the length of geodesies on a Riemann surface and the eigenvalues of the 
Laplacian operator defined on it This duality is encapsulated in the so called Selberg 
trace formula, which has a strong similarity with the Riemann explicit formula relating the 
zeros and the prime numbers. The Riemann zeros would correspond to the eigenvalues, 
and the primes to the geodesies. This classical versus quantum version of the primes and 
the zeros was also at the heart of the so called Quantum Chaos approach to the RH (see 
later). Quite independently of the Selberg work, Montgomery showed that the Riemann 
zeros are distributed randomly and obeying locally the statistical law of the Random Matrix 
Models (RMM) Q. The RMM were originally proposed to explain the chaotic behaviour 
of the spectra of nuclei but it has applications in another branches of Physics, specially 



in Condensed Matter 



13] . There are several universality classes of random matrices, and 



it turns out that the one related to the Riemann zeros is the gaussian unitary ensemble 
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(GUE) associated to random hermitean matrices. Montgomery analytical results found an 



impressive numerical confirmation in the works of Odlyzko in the 80 's, so that the GU 



as applied to the Riemann zeros is nowadays called the Montgomery-Odlyzko law 14]. 



law, 



It is worth to mention that the prime numbers, unlike the Riemann zeros, are distributed 
almost at random over the set of integers. Indeed, it is believed that one can find arbitrary 
pairs of nearby odd numbers n, n + 2, as well as pairs arbitrarily separated. The only thing 
known about the distribution of the primes is the Gauss law according to which the n th prime 

n 

p n behaves asymptotically as p n ~ n/\ogn [1]. This statement is called the Prime Number 
Theorem proved by Hadamard and de la Valle-Poussin in 1896. If the RH is true then the 
deviation from the Gauss law is of order y/n logn. The analogue of the Gauss law for the 
imaginary part of the Riemann zeros (called it E) is given by the Riemann law where the 
n th -zero behaves as E n ~ 2im/ logn. Hence, large prime numbers are progressively scarced, 
while large Riemann zeros abound. 

An important hint suggested by the Montgomery-Odlyzko law is that the Polya-Hilbert 
Hamiltonian H must break the time reversal symmetry. The reason being that the GUE 
statistics describes random Hamiltonians where this symmetry is broken. A simple example 
is provided by materials with impurities subject to an external magnetic field, as in the 
Quantum Hall effect. 

A further step in the Polya-Hilbert-Montgomery-Odlyzko pathway was taken by Berry 
, |4J . who noticed a similarity between the formula yielding the fluctuations of the number 
of zeros, around its average position E n ~ 27m/ logn, and a formula giving the fluctuations 
of the energy levels of a Hamiltonian obtained by the quantization of a classical chaotic 
system 15|. The comparison between these two formulas suggests that the prime numbers 



p correspond to the isolated periodic orbits whose period is log p. In the Quantum Chaos 
scenario the prime numbers appear as classical objects, while the Riemann zeros are quantal. 
This classical/quantum interpretation of the primes/zeros is certainly reminiscent of the one 
underlying the Selberg trace formula mentioned earlier. One success of the Quantum Chaos 
approach is that it explains the deviations from the GUE law of the zeros found numerically 
by Odlyzko. The similarity between the fluctuation formulas described above, while rather 
appealing, has a serious drawback observed by Connes which has to do with an overall sign 
difference between them jl6]. It is as if the periodic orbits were missing in the underlying 
classical chaotic dynamics, a fact that is difficult to understand physically. This and other 
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observations lead Connes to propose a rather abstract approach to the RH based on discrete 



mathematical objects known as adeles 



161 ] . The final outcome of Connes work is a trace 



formula whose proof, not yet found, amounts to that of a generalized version of the RH. 
In Connes approach there is an operator, which plays the role of the Hamiltonian, whose 
spectrum is a continuum with missing spectral lines corresponding to the Riemann zeros. 
We are thus confronted with two possible physical realizations of the Riemann zeros, either 
as point like spectra or as missing spectra in a continuum. Later on we shall see that 
both pictures can be reconciled in a QM model having a discrete spectra embedded in a 
continuum. 



II. SEMICLASSICAL APPROACH 



In 1999 Berry and Keating on one hand 171. Il8j|. and Connes on the other 16|, proposed 



that the classical Hamiltonian H = xp, where x and p are the position and momenta of a ID 
particle, is closely related to the Riemann zeros. The classical trayectories of the particle are 
hyperbolas in the phase space (x,p), hence one should not expect a discrete spectrum even 
at the semiclassical level (see fig. [TJ To overcome this difficulty, Berry and Keating proposed 
to restrict the phase space to those points where \x\ > l x and \p\ > l p , with l x l p = l-nh. The 
number of semiclassical states, J\f(E), with energy between and E is given by the allowed 
area in phase space divided by h = 2n {h = 1) (see fig. [I]) 

Af(E) = ^-(log^--l)+l (1) 
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Eq.([T]) coincides asympotically with the smooth part of the formula that gives the number 
of Riemann zeros in the same interval. This result is really striking given the simplicity 
of the Hamiltonian and the asumptions made. On the other hand Connes started on the 
same classical Hamiltonian H = xp but constrained the phase space to those trayectories 
satisfying \x\ < A, \p\ < A, with A a cutoff which is sent to infinity at the end of the 
calculation. The number of semiclassical states is given by 

*f{E) = *L log A - ^ flog |; -l) (2) 

The first term describes a continuum of states in the limit A — > oo, while the second term 
coincides with minus the average position of the Riemann zeros (TTT) . This result lead Connes 




FIG. 1: la) a classical trayectory of the Hamiltonian H = xp. The regions in shadow are the 
allowed phase space of the semiclassical regularizations of H = xp considered by: lb) Berry and 
Keating, lc) Connes and Id) Sierra. 

to propose the missing spectral interpretation of the Riemann zeros described earlier. 



2fJ, is that 



A third possible regularization of the xp model, proposed in references 
l x < x < A, which leads to the following counting of semiclassical states, 

This result agrees with the asymptotic part of (jSJ), meaning that there is a continuum 
spectrum but the possible connection with the Riemann zeros is lost. The main advantage 
of the latter regularization is that it arises from a consistent quantization of H = xp unlike 
the two previous semiclassical regularizations. 

III. QUANTIZATION OF xp AND l/(xp) 

The classical Hamiltonian H = xp can be consistently quantized in two cases de pen ding 



on the choice of the domain in the x coordinate: 1) < x < oo and 2) l x < x < A 



20. 



2ii 
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In the first case H is essentially self- adjoint, while in the second it admits a one parameter 
self-adjoint extension. We shall consider the latter case. To do so one first define the normal 
ordered operator [l^] 

#o = ^(xp + px) = + ^) (4) 

where p = —id/dx. The formal eigenfunctions of (j3J) are 

Ms) = -^=iEi 1 < x < N (5) 

where we have normalized l x = 1 and A = N. One can show that (J4]) is self-adjoint if the 
wave functions satisfy the boundary condition 

e^W = iV 1/2 MN) (6) 
where the angle 9 parameterizes the self-adjoint extension of H . Imposing ([6]) on (jSJ) yields 

N iE = e * e ^ 

which determines the eigenvalues of H 

2tt / 9 \ 

£ "=toilv(" + ^j' - = ".±1.- W 
This equation agrees with the semiclassiclassical result (j3]). In the particular case where 
9 = 7r, the spectrum (jSj) becomes symmetric around 0. Notice that the zero eigenvalue is 
excluded. Another way to derive this result is by considering the inverse of the operator (j4j). 
This can be done as follows, 

Ho = ~(xp + px) = £2 p £2 — » H^ 1 = x _ 2 p~ x x~2 (9) 

where x 1 ^ 2 is well defined since x > 1. The inverse of the momenta operator p is given by 
the ID Green function 

p -i = » sign(s-s') ^ (1Q) 
2 \Jxx' 

where sign(x — x') is the sign function. The Schroedinger equation associated to (jU]) is 
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i f N ^ sign(x x') ^ = 1 
* Ji Vxx' & 



which in terms of the new wave function 



becomes 



<f> E (x) = x 2 ^e(x) (12) 



iE f N 

x 4>e(x) / dx sign(x — x')4>e(x) = (13) 

2 J\ 



whose solution is 



A^ = -l, 4>{x) = -£- (14) 



C 



x 



Hence we recover eqs.flS]) and (jTJ) in the particular case where 9 = ir. On the other hand, 
eq. (fl3|) looks as the eigenvalue equation of yet another Hamiltonian that we shall discuss 
next. 



IV. RELATION WITH THE RUSSIAN DOLL BCS MODEL 



In reference [22l. l23l. I24J] it was defined an extension of the BCS model of superconductivity, 
called the Russian Doll model, whose Hamiltonian, when restricted to the one body case 
becomes 

H RD (x, x') = e(x) 8(x — x') — - (g + i h sign(x — x')) (15) 

where e(x) represents the energies of pairs of electrons occupying time reversed states in the 
band 1 < x < N, g is the standard BCS coupling constant and h a coupling that breaks the 
time reversal symmetry. The eigentates and eigenfunctions of ( JTBl are given by 

N-E RD \ lh g + ih C 



1 - Erd J g -ih (x — E 



RD 



Comparing (fl6l) with pijl . one obtains the following map between the eigenstates of the xp 
model and the RD model 
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E RD = <-> E ± (17) 
h*->E 
£ «-> tan(0/2) 

«-> X~ l ' 2 ll)E 

in the particular case where g = and then 6 = n. One can add a g coupling term in the 
definition (jHJ) of the operator if - , in which case the correspondence between the RD model 
and the xp model will cover all the self-adjoint extensions of xp. The RD model provides an 
example where the renormalization group, instead of ending at fixed points, run in cycles 



2a 



26, 



271 ] . In the RD model, the coupling that runs periodically under the RG is g, with a 



period equal to 2ir/h, while the coupling h remains invariant. This fact in turn implies the 
existence of several bound states whose number is given by n = 2n/hlogN. If we replace 
h by E, the latter equation becomes n = 2ir / 'ElogN which coincides with the number of 
eigentates of the xp model. 

Incidentally, we would like to mention that the field theory realizations of the cyclic 

3o| . are at the origin of LeClair's approach 



Renormalization Group, of references 



to the RH 



28, 



29 



3l| . In this reference the zeta function on the critical strip is related to the 



quantum statistical mechanics of non-relativistic, interacting fermionic gases in Id with a 
quasi-periodic two-body potential. This quasi-periodicity is reminiscent of the zero temper- 
ature cyclic RG of the quantum mechanical Hamiltonian of [^(J , but the general framework 
of both works is different. 



V. H = xp WITH INTERACTIONS: GENERAL MODEL 

The previous results establish an interesting correspondence between two apparently dif- 
ferent models which also suggests a way to add interactions to the xp model. Indeed, the 
interacting term of the RD Hamiltonian (fl"5l) that is proportional to the coupling constant 
g is basically a proyector operator \BCS){BCS\, with wave function (x\BCS) = 1 Vx. As 
we said above, adding that term to the inverse Hamiltonian ([9]) would give rise to a 9 ^ n 
term associated to the self-adjoint extensions of xp. Instead we want to add an interacting 
term that reflects the existence of two boundaries in the BK regularization of the xp model. 
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The simplest possibility is to define 



H 



-i 



H^ + -(\^ a )(^ b \-\^ b )M) 



(18) 



where il) a ,b are two wave functions whose properties will be specified below. The matrix 
elements of ([TBI read 



tt-\ ( n i sign(x - x') + a{x)b(x') - a{x')b(x) , 

H [X, X ) — ; , 1 < X, X < JM 



V. 



XX 



(19) 



where 



a(x) b{x) 
x&> MX) = x^ 
r-1 



(20) 



are real functions, which guarantee that H~ is a hermitean and antisymmetric matrix, so 
that its eigenvalues are pairs of real numbers E, —E ^ 0. A simplified version of (fT9l) is 
obtained by choosing b(x) = 1. The latter model will be denoted as type I, while the former 
as type II. For these models to be well defined in the limit iV — > oo we impose the following 
normalization conditions 



f(x) n 

dx LL-L\ < oo, (n= 1,2), 



x 



(21) 



a type I 
a, b type II 

A nice feature of the Hamiltonian (fT9j) is that the Scroedinger equation is exactly solvable, 
the reason being that it is equivalent to a first order differential equation supplemented with 
a boundary condition. We shall next present the results obtained in reference 20(. First of 
all, the eigenenergies E of (fUl satisfy the equation 



T N {E)+T N {-E) = 



(22) 



where !Fn(E) is a Jost like function whose expression will given below in the limit iV — > oo. 
In that limit the eigenf unctions of the type II model are given by 



X 



1 

1/2-iE 



Coo + 



dyy 



-iE 



db 



( da ^ 
\dy dy 



A 



(23) 
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where A, B, are integration constants that depend on E. In the limit x — > oo the 
functions a(x),b(x) vanish sufficiently fast so that the wave function (123]) is dominated by 
the first term, i.e. 



C 

x»i K ' x l/2-*E 



(24) 



It turns out that is given by the Jost function (up to a constant which can be taken as 
1), 

C OQ {E)=J r {E) (25) 

Hencefore the energies where J~(E) is non zero correspond to delocalized states which behave 
asymptotically as the eigenf unctions of the unperturbed Hamiltonian H Q = (xp+px)/2. For 
these states the ratio J r (E)/J-'(—E) gives the scattering phase shift. On the other hand, 
vanishes whenever ^{E) does. These energies correspond to localized states with a finite 
norm. In summary, the spectrum of H consists of a continuum formed by those energies 
where T{E) ^ 0, plus a discrete part given by the real zeros of T{E) (see fig. [2]). 

Moreover, using the hermiticity of H one can show from eq. (!2ip that J- co {E') does not 
have zeros with lm(E) > 0. 

T{E) = =>• Im E < (26) 

The real zeros of J~{E) correspond, as explained above, to localized states, while the complex 
zeros below the real axis correspond to resonances. These results are summarized in table 
1. 



Eigenstate 


Coo 


HE) 


Eigencondition 


Delocalized 
Localized 


= 


= 


ATiE _ HE) 
JV ~ T(E)* 

F{E) = 



Table 1.- Classification of eigenstates of the model. 

Before giving the expression of T{E) for the type I and type II model we shall introduce 
some definitions. First of all let us define the Mellin transform of a (similarly for b(x)). 



a(t) 



dx x 



-l+it 



a[x) 



a(x) 



dt 
2^ 



x u a(t) 



(27) 
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Bound states 



Scattering states 




FIG. 2: Pictorial representation of the spectrum of the model. The bound states are the points 
where J~(E) = 0, which are embedded in a continuum of scattering states. 

The reality of a(x) implies 

a*{t) = a(-t), t e R (28) 
Condition n — 1 in eq. (l2Tl) amounts to 



|a(0)| < oo (29) 
while condition n = 2 in eq. (l2Tj) is equivalent to 

r dx = r *l \z( t )\ 2 < oo, (30) 

Ji % y_ oc 2tt 

The function a(t) is in fact the Fourier transform of a(x) in the variable q = log x, which takes 
values in the interval (0, oo). In terms of q, a(x) is a square normalizable and causal function, 
in which case a(t) has interesting analytic properties by a theorem due to Titchmarsh 32]. 
This theorem states that under the previous conditions a(t) is analytic in the complex 
upper-half plane and satisfies the formula 



r , , dt a(t) , 
a(z) =P / — zen 31 

tlT t — Z 



oo 



where P denotes the Cauchy principal value of the integral. To prove this formula one 
uses the fact that a(t) has no singularities in the upper-half plane and that the contour on 
integration on the circle \z\ = R, Im z > vanishes since limM._ >00 | a(z) | = 0. For later 
purposes let us define the new function 

a(t) = l ±a(t), (32) 
and similarly h(t), whose properties follow from those of a(t) namely: 
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Reality: 



a*(t) = a(-t), t e R (33) 



Regularity: 



lim if^ll < oo =^ | a ( ) | = (34) 



t-»o t 



Normalizability 



dt \a(t)\ 2 

1 < oo, (35) 



-co 2tt t 2 



Analiticity 



a(z) = P —-AJ--P zen (36) 



oo in t - z > /_ 00 in t 



This eq. implies that a(z) is an analytic function in the upper-half plane which con- 
verges towards a constant value in a circle of infinite radius given by the last term of 

(EED. 



The next definition involves the product of two analytic functions f(t) and g(i) in the 
upper half-plane: 

(f *g)(z) = f(z) g(-*) + r - i{t \ Si ~ Z \ z E R (37) 

where the integration is understood in the Cauchy sense as in eqs. (I31j) and fl36l) . One can 
show that (f -k g)(z) is an analytic function in the upper half-plane provided f(z) g{—z) is 
well behaved, which seems to be the case in all the examples we have analized. The analytic 



extension of (f * g) (z) to the lower half-p 
this product we shall define the function 



ane will have in general singularities. In terms of 



33|: 



StM = {f* S ){z)-{f*s){0) (38) 
which satisfies the following conditions 
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Reality: if f and g verify (|33|) then 



S{, s {z) = 5 f , g (-z), zeiR (39) 

• Regularity 

5 fig (0) = 1 (40) 

• Shuffle relation 

StM + Sui-z) = 2 f(z) g(-s) (41) 

The notation "shuffle" is borrowed from the Theory of multiple zeta functions, as 
explained later on. 

After these definitions we can finally give the expression of the Jost function J-{E) in terms 
of the potentials a(i) and b(i). For the type I model it reads 

F{t) = 1 + 2 a(t) + S a>a (t) (42) 

while for the type II model it is 

F(t) = 1 - 5 a,b(t) + S h}lL (t) + S^(t)S hih (t) - S &}h (t)S h}lL (t) (43) 
From the properties of the S"-functions one can easily derived: 

• Reality: 

F*(t) = T{-t) z G 1R (44) 

This condition guarantees that the ratio J r (E)/J-'(—E) appearing in the eigenvalue 
eq. (1221) is indeed a phase. 

• Regularity 



^(0) = 1 
13 



(45) 



This condition implies that the zero energy is not an eigenvalue of the Hamiltonian 
H, which was the assumptions made by defining it in terms of its inverse. 

Let us next consider the two models separately. 
A. Type I model 

The Jost function fj42l can be expressed as 

f(t) = l + 2a(f) + (a*a)(J) (46) 
where we used that a(0) = 0, which implies 

(a*a)(0) = 0. (47) 
An important consequence of eq. (l4~6l) is the positivity of the real part of ^(t), 

Re F(t) = |1 + a(t)| 2 > (48) 

which imposes a strong contraint on the functions allowing a QM interpretation as Jost 
functions of the type I model. In particular eq. (j48|) excludes the zeta function ((a — it) for 
1/2 < a < 1, but not the case where a > 1 as we shall see later on. 

The ^-product defined in eq. (|37l) is non commutative and non associative. Nevertheless 
it behaves nicely respect to the identity function 1, namely 

1*1 = 1 (49) 

and 

a*l + l*a = 2a (50) 

where a is an analytic function in the upper half-plane satisfying eq. (1361) . Using these 
two equations, the Jost function ( H6|) can be expressed as the square, with respect to the 
^-product, of a single function, i.e. 

jr = f*f ; f=l + a (51) 
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so that f is the ^-square root, of T . Using eq. (l5T]) one can easily prove that T{f) does not 
have zeros in the upper half-plane. Indeed, write f * f as 



(f*f)(s)= / % m t ^\ ImoO (52) 



— oo 



then 

Re(f*0W=£;fj^L>0, , = , + * (53) 

A simple but illustrative example of the theory developped so far is provided by the step 
potential, 



.Oi, 1 < X < X\ 

a(x) = { (54) 
0, x\ < x < oo 

which yields 

a=|(xf-l), a *a=|(l-*f), (55) 

and the Jost function 

m = ! + a -l^pl {x 'f - l). (56) 

Fig. [3] shows an Argand plane representation of the real and imaginary parts of f)56p for 
several values of a±. For a± = 1 the function !F(t) vanishes at the values 

(2n+l)vr . , 

E n = y — J -, n = 0,±l,... (57) 

log^i 

describing an infinite number of bound states. The remaining values of E correspond to 
delocalized states. All the zeros of T{E) lie on the real axis for a x = 1 and below it for 
at ^ I. 



B. Perturbative solution of the type I model 

The next problem we address is: given a function J-{t), satisfying the analyticity, reality, 
regularity and positivity conditions described above, which is the potential, or potentials, 
a(t), verifying eq. (|46|) ? In this paragraph we shall give a perturbative method to construct 
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FIG. 3: Plot of the real and imaginary parts of J-{E), as given by eq. (|56p . for the choices a\ = 
0.2, 0.4, 0.6, 1. At a% = 1 the circle passes through the origin. 

one of those potential in terms of a series which converges under certain conditions placed 
on the function J-(t). Let us first make the change a — > —2a in eq. (l46|) which becomes 

1 - T 

a = g + a*a, g = - (58) 
Iterating (1551) generates the series expansion 

a = g + g*g + g*(g*g) + (g*g)*g + ... (59) 

At order g™ one gets all admissible bracketings, whose number is given by C n _i, where C n 
is the Catalan number 



n + l \ n I V^r n A l l 



whose asymptotic behaviour is also described. To investigate the conditions for convergence 
of (159|) . let us suppose that g is given by the absolute convergent series 

oo oo 

g(*) = X>< |g(*)|<£)|</n|<00 (61) 
n=l n=l 

Using (1371) one finds 
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(g * g) (*) = J2 9n + 2 E 9n 9m (n/mf (62) 

n=l n>m 

so that 

oo 

|(g*g)(*)| < £l^n|) 2 , VtGR (63) 



n=l 



Similarly one finds that 



|a(f)|< j^CU j^bml) , (64) 



n=l \m=l 

which, according to (I60I) . converges provided 



m=l 



This condition is sufficient for convergence of the series fl59l) but it is not necessary. Eq. fl65|) 
implies 



\T{t) - 1| < 1, Vt G R (66) 

However the converse is not true (i.e. eq. fl66l) does not imply fl65l)). We believe that fl66|) 
also guarantees the convergence of (!59|) . but this guess needs to be proved. Observe also 
that fl66l) implies that Re .F(i) > 0, which is also a necessary condition for the existence of 
the potential a(t). An application of these results is: 



C. The zeta function with a > 1 

For a > 1 let us consider the function 

F(t)=C((a-it), C=-?- (67) 

where the constant C guarantees the normalization condition .F(O) = 1. The values of the 
constants g n appearing in flBTl) are given by 

1 — C C 
9i = ^ r , 9n>i = - I ^ (68) 
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and the convergence criteria (I65I) yields 



Ew = — <r^- 2 (69) 

m=l 

where the latter conditions is satisfied if 



a > a c = 1.72865, C(0 = 2 (70) 
Moreover condition (1661) can be checked numerically i.e. 

| C( CT~^ -1|<1, V(T>1 ^ Wg]r 

so we expect that the series (I59"j) will also converge for any value a > 1. Fig. H] displays the 
real and imaginary parts of a(t) for the case a = 2 obtained by the sum of the first terms of 
eq. fl59|) . The convergence towards a finite value is clear. 

The series (l59l) contains also an interesting analytical structure, which can be seen from 
the star product of two zeta functions, 

oo oo 1 

C(" - «) = E ^ - (C * 0(*) = C(2a) + 2 £ — (72) 

n=l n>m>l 

The double sum series of the RHS is equal to a Euler-Zagier zeta function for two variables 
which is a generalization of the zeta function. Multivariable versions of this function have 
attracted much attention in various fields, as knot theory, perturbative quantum field theory, 



etc (see |34j, [35( and references therein). 



D. The zeta function at a = 1/2 

The results obtained so far suggests that the Riemann zeta function on the critical line 
C(l/2 — iE) could perhaps be realized as the Jost function of the model. This idea is 
motivated by the scattering approa ch p ionered by Faddeev and Pavlov in 1975, and has 
been followed by many authors 36|, |37J, |38|, |39j . An important result is that the phase of 



((1+it) is related to the scattering phase shift of a particle moving on a surface with constant 
negative curvature. The chaotic nature of that phase is a well known feature. Along this line 
of thoughts, Bhaduri, Khare and Law (BKL) maded in 1994 an analogy between resonant 
quantum scattering amplitudes and the Argand diagram of the zeta function ((1/2 — it), 
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FIG. 4: Real and imaginary parts of the potential a(i) for a = 2 in the interval t E (10, 25). 



where the real part of £ (along the x-axis) is plotted against the imaginary part (y-axis) [40 ]. 
The diagram consists of an infinite series of closed loops passing through the origin every 
time C(l/2 ~ it) vanishes (see fig. [5]). This loop structure is similar to the Argand plots 
of partial wave amplitudes of some physical models with the two axis being interchanged. 
However the analogy is flawed since the real part of ({1/2 — it) is negative in small regions 
of t, a circumstance which never occurs in those physical systems. 

In fact, the loop structure of the models proposed by BKL is identical, up to a scale 
factor of 2, to the model with the potential (|54|) (see fig. ([3D), where the loops representing 
F(E), for d = 1, are circles of radius 1/2, centered at x — 1/2. For general models of 
type I, the loops are not circles but the real part of J-'i(E) is always positive (see eq. (l4"8l) ). 
and therefore they can never represent ^(1/2 — iE). Incidentally, this constraint does not 
apply to the models of type II, where Re F(E) may become negative. This suggests that 
((1/2 — iE) could indeed be the Jost function J-'(E) of a type II model for a particular choice 
of a and b. If this were the case then the Riemann zeros would become eigenenergies of the 
Hamiltonian realizing in that manner the Polya and Hilbert conjecture which may also give 
hints into the solution of the Riemann hypothesis. A complete answer to this problem is 
not yet known but we shall present below some encouraging results along this direction. 

The first step is to recover quantum mechanically the smooth approximation to the 
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Riemann zeros. This approximation is equivalent to the following condition 



1 + = (73) 

where 

e mE) = K^rM^ ( 74 ) 

1 U 2 ) 

The function 9(E) gives the phase of the zeta function, e.g. 

C(l/2 - it) = e im (75) 

while is the Riemann- Siegel zeta function which is real and even due to the duality 
relation satisfied by the zeta function (e.g. C(V^ — it) = e 2ld ^ <^(l/2 + it)). The reason 
that ( 1731) is a good approximation to the location of the zeros, can be seen in fig. [5] which 
plots the real and imaginary parts of C(l/2 — iE). Observe that in the vicinity of a zero, 
the curves cut the imaginary axis where 9(E n ) = ir(2n — l)/2 so that cos 9(E) = 0, which is 
nothing but eq. (I73l) . The value of E n satisfies 

n = —— + - (76) 

7T Z 

which asymptotically coincides with eq. jT| up to a constant term. This result can be 
obtained choosing the following potential in the type I model 

a(x) = -2 Sin(2 " x) (77) 

<x 



whose Mellin transform (1271) gives 

a( t ) = -je 2i9 ^ +0(l/t) (78) 
up to 1/t terms. The corresponding Jost function is 

F(E) = 2(1 + e 2i9{E) ) + 0(1/ 'E) (79) 

whose zeros agrees with (1731) asymptotically. Hence the smooth approximation to the Rie- 
mann zeros can be obtained asymptotically by the potential (1771) . Why this potential is able 
to yield this result? One may suspect that it must implement at the quantum level the BK 
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FIG. 5: Left: real and imaginary parts of with s = 1/2 — iE and E € (0,50). Right: same as 
before for Ch(s) = (s — l)((s)/s. 

boundary conditions of the semiclassical approach. Indeed, let us show how this works in 
detail. The potential (ITT)) corresponds to the wave function 

, , . „ sin(27rx) 

i) a {x) = -2 —± 1 80 

x 

which satisfies the equation 

H*il> a (x)= ((2nx) 2 -^j^ a (x), (81) 

where H = -Jxp^fx is the Hamiltonian (jUJ). Dropping ip a in both sides and replacing H by 
xp, one obtains a classical version of ( JSTT) . 



{xpf = (2nxf -\^ V = ±2Jl- (82) 

which describes a curve in phase space that approaches asymptotically the lines p = ±2ir. 
We shall identify these asymptotes with the BK boundary in the momenta \p\ = l p . Recall 
on the other hand the boundary condition x > l x = 1, which combined with the previous 
identification reproduces the Planck cell quantization condition, 

l p = 27r, l x — 1 ==}- l p l x = 2n (83) 
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The BK condition \x\ > l x has already built in into the model and this is also reflected in 
the state ipb(%) — Vv 7 ^ which is concentrated near the position l x = 1. In the model we 
have proposed, the two BK boundary conditions are realized asymmetrically, as opposed to 
the semiclassical model. It would be desirable to have a more symmetric treament of them. 
This indeed can be done and the results will be presented elsewhere. 

The next natural step is to see wether the zeta function C(l/2 — iE) can be realized as 
the Jost function of the model. The analyticity properties of the Jost functions implies that 
T{E) must be of the form 

F{E) = C^-^- C(l/2 - iE) (84) 

which does not have poles in the upper-half plane (C and fi > are normalization constants). 
In fig. [5]we plot an example with \x = 1/2. Since the real part of ({1/2 — iE), and thus 
of ^F(E), become negative in small regions of E one is forced to consider the type II model 
with two non trivial functions ip a and ipb- The problem of finding ip a ^ is rather non trivial. 
One can in principle fix one of them, say ipb, and try to solve for ip a as a function of the 
Jost function (I84p and ipb- in the case of the zeta function ((a — iE), with a > 1, we 
were able to solve this problem perturbatively thanks to the fact that the zeta function is 
bounded. However for a = 1/2 the zeta function is unbounded which may lead to problems 
of convergence. In any case, it seems clear that one needs further physical insights to 
make progress into this difficult problem. As we said above, one needs a more symmetric 
treatment of the coordinates x and p, and a clearer physical interpretation of the wave 
functions ip a ,b- Another important ingredient to be implemented is the duality symmetry of 
the zeta function which in the present formulation of the model is not manifest but which 
is expected to play a central role. 

In summary, we have presented in this work an interacting version of the xp Hamiltonian 
which may ultimately lead to a spectral realization of the Riemann zeros, as suggested long 
ago by Polya and Hilbert. The main ingredient of the model is the non local character of 
the interaction in terms of two potentials which are the quantum version of the semiclassical 
phase space constraints of Berry and Keating. The generic spectrum of the model consists 
of a continuum of eigenstates in the thermodynamic limit which may also contain bound 
states embedded in it. We conjecture the existence of potentials giving rise to the Riemann 
zeros as the discrete spectrum embedded in the continuum. If this were the case that would 
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resolve the emission versus absortion spectral interpretation of the Riemann zeros. This 
would also open the way to a better understanding of the Riemann hypothesis. We have 
also pointed out the need to implement in an explicit way the duality properties of the zeta 
function, which implies a more symmetric treatment of the x and p coordinates as in the 
semiclassical model. 
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